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Abstract 

The Mars flyby of C/2013 Al (Siding Spring) represented a unique opportu¬ 
nity for imaging a long-period comet and resolving its nucleus and rotation 
period. Because of the small encounter distance and the high relative ve¬ 
locity, the goal of successfully observing C/2013 Al from the Mars orbiting 
spacecrafts posed strict accuracy requirements on the comet’s ephemerides. 
These requirements were hard to meet, as comets are known for being highly 
unpredictable: astrometric observations can be significantly biased and non- 
gravitational perturbations affect comet trajectories. Therefore, even prior 
to the encounter, we remeasured a couple of hundred astrometric images ob¬ 
tained with ground-based and Earth-orbiting telescopes. We also observed 
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the comet with the Mars Reconnaissance Orbiter’s High Resolution Imaging 
Science Experiment (HiRISE) camera on 2014 October 7. In particular, these 
HiRISE observations were decisive in securing the trajectory and revealed 
that out-of-plane nongravitational perturbations were larger than previously 
assumed. Though the resulting ephemeris predictions for the Mars encounter 
allowed observations of the comet from the Mars orbiting spacecrafts, post¬ 
encounter observations show a discrepancy with the pre-encounter trajectory. 
We reconcile this discrepancy by employing the Rotating Jet Model, which 
is a higher hdelity model for nongravitational perturbations and provides an 
estimate of C/2013 Al’s spin pole. 

Keywords: Comets, Comets, dynamics. Data reduction techniques. Orbit 
determination 


1. Introduction 


On 2014 October 19 long-period comet C/2013 Al (Siding Spring) ex¬ 
perienced an exceptionally close encounter with Mars at 140496.6 ± 4.0 km 
and a relative velocity of 55.963249 ± 0.000028 km/s (la formal uncertain¬ 
ties). While an impact between the nucleus of C/2013 Al and Mars had been 


ruled out by earlier observational data, early predictions by Vaubaillon et al 


(2014) and Moorhead et al. (2014) suggested that the dust in the comet’s tail 


could have posed a signihcant hazard to the spacecrafts orbiting Mars. Later 


studies (Farnocchia et ah, 2014 Kelley et ah, 2014 Tricarico et ah, 2014; Ye 


and Hui, 2014) used observations of C/2013 Al to estimate the dust produc- 
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tion rate and concluded that spacecrafts were unlikely to be affected by the 
Mars flyby of C/2013 Al. Moreover, these studies predicted the fluence peak 
timing as well as the direction of the incoming dust to inform the rephasing 
of the spacecraft orbits. As expected, no damage was reported. 

Along with the potential hazard, the close encounter between Mars and 
C/2013 Al represented a unique scientihc opportunity: never had a long- 
period comet been observed at such small distance. The instruments on 
board of the Mars orbiting spacecrafts could be used to collect valuable 
observations of the comet and for the hrst time allow one to resolve the 
nucleus and rotation period of a long-period comet. 

The goal of observing C/2013 Al during the Mars encounter posed strict 
requirements on the ephemeris accuracy. While the Mars Reconnaissance Or- 
biter (MRO) High Resolution Imaging Science Experiment (HiRISE) camera 


has a very large held of view capability (McEwen et ah, 2007), only a small 


portion of its held of view, 4 mrad x 4 mrad, could be used for comet ob¬ 
servations because of data handling limitations. This translates to a 280 km 
error tolerance at a distance of 140,000 km. Meeting these requirements was 
not an easy task as comets are famously difficult to predict]^ As a matter 
of fact, because of comet tails and outgassing, comet astrometric observa¬ 
tional data are not as clean as for other bodies, e.g., asteroids. Moreover, 
nongravitational perturbations often limit the capability of providing accu- 


^ Comets are like cats; they have tails, and they do precisely what they want, David H. 
Levy 
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rate ephemerides (Yeomans et al., 2004): these perturbations can be hard 
to model and estimate due to the lack of information and yet signihcantly 
affect comet trajectories. 

Comet ephemeris challenges associated with the Deep Space 1 mission’s 


encounter with Comet 81P/Borrelly in 2001 (Rayman, 2002) led to the de¬ 


velopment of a rotating jet model for nongravitational accelerations (Chesley 


and Yeomans, 2005), which provided far superior comet ephemeris predic¬ 


tions, albeit ex post facto. Later, the ephemeris development efforts for 
subsequent NASA comet flyby missions, such as Deep Impact, Stardust, 
EPOXI and Stardust-NExT, went smoothly, in large part because of the im¬ 


provements in nongravitational acceleration modeling (see, e.g., Chesley and 


Yeomans, 2006). 


As we shall see, the technical approaches developed during these hve 
comet flyby missions were invaluable in deriving an accurate ephemeris for 
the Siding Spring encounter of Mars. A persistent challenge in such cases 
revolves around understanding how astrometric bias affects the measured 
nucleus position due to the effects of an asymmetric coma. Complicating 
the situation is the possible presence of nongravitational accelerations that 
could be affecting the trajectory in ways that are difficult to distinguish from 
astrometric bias. Separating these two sources of error, nongravitational ac¬ 
celerations and astrometric bias can be helped by carefully examining fits 
and predictions associated with subsets of the astrometric data set and by 
taking special care in the astrometric reduction to measure the peak of the 
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condensation in the target’s point spread function. Fortunately, as the time 
to encounter shortens the problem tends to resolve itself due to the combi¬ 
nation of longer data arcs, which tends to improve nongravitational acceler¬ 
ation modeling, and shorter mapping times, which reduce the negative effect 
of modeling errors. 


2. Analysis prior to the encounter 


At the end of March 2014, C/2013 Al went into solar conjunction thus 
becoming unobservable from Earth. The trajectory obtained with the pre¬ 
conjunction observed arc was JPL solution 46, which is fully described in 


Farnocchia et al. (2014). At the end of May 2014 C/2013 Al emerged from 


the Sun, thus leading to additional astrometric observations and correspond¬ 
ing orbit updates. However, these orbit updates showed a quite erratic be¬ 
havior. For different orbital solutions. Fig. shows the predicted locations 
on the 2014 October 19 encounter 6-plane, which is the plane containing the 
Mars center at the coordinate origin and normal to the asymptotic incoming 


velocity of the comet Voo (Kizner, 1961). The projection of the heliocentric 


velocity of Mars on the 6-plane is oriented as —while the ^ coordinate 
completes the reference frame, i.e., ^ = Voo x C/|voo| (for more details see 


Valsecchi et al., 2003). The sequence of orbit updates from solution 46 to so¬ 


lution 56 was statistically inconsistent and raised questions about the ability 
to provide sufficiently accurate encounter predictions. 
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Figure 1: Mapping of the C/2013 A1 Itr position uncertainty on the 2014 October 19 
6-plane for sequential orbital solutions. Dashed ellipses are for earlier gravity-only solu¬ 
tions computed before the astrometry selection and remeasurement process described in 
Sec. |2.1| solution 46 includes data until 2014 March 31, while solution 56 includes data 
until 2014 August 11. Solution 95 and 97 are the two solutions computed before the 2014 
October 7 HiRISE astrometry, while solution 101 is the final ephemeris delivery prior to 
the encounter. Finally, solution 102 is computed employing the Rotating Jet Model and 
using the full data arc through end of April 2015. 
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2.1. Astrometry 

One of the sources of the erratic behavior shown in Fig. [T] is the poor 
quality of the astrometry used to compute the orbital solutions. Astrometry 
of active comets is often affected by an asymmetric distribution of light from 
the coma surrounding the nucleus. If a model designed to £t the symmetric 
light distribution from a point source is used on an active comet, the fitted 
photocenter can be dependent on the size of the synthetic aperture utilized. 

This effect was quite obvious in the astrometry of comet P/Wild 2, ob¬ 


tained in support of the Stardust mission’s flyby of that comet (Tholen and 


Chesley, 2004). As the size of the synthetic aperture increased, the measured 


photocenter moved in the tailward direction. At least in the case of P/Wild 
2 and over the range of synthetic apertures tested, the offset showed a fairly 
linear dependence on synthetic aperture size, so the position of the cometary 
nucleus was computed by linearly extrapolating the fitted coordinates to a 
theoretical zero-sized synthetic aperture. The resulting astrometry enabled 
the uncertainty in the position of P/Wild 2 to be shrunk by a factor of three 
prior to the spacecraft encounter. 

The same tailward bias was observed in the astrometry of C/2013 Al 
as the size of the synthetic aperture increased, therefore one should take 
the astrometric position corresponding to the extrapolation of the trend to a 
theoretical zero-sized aperture or at least to the brightness peak. However, we 
had no control nor information on most of the astrometry that was reported 
to the Minor Planet Center, some of which obtained by amateur observers. 
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To obtain more reliable orbital solutions, we went through a strict pro¬ 
cess of astrometry selection and remeasurement. First, we selected the few 
astrometric observations obtained by professional surveys: Pan-STARRS 1, 
Kitt Peak, and Siding Spring Survey. It is worth pointing out that most of 
the time C/2013 Al was at high negative declination, and therefore was not 
observable from the Northern Hemisphere where most professional surveys 
are located. Then, we added observations that we obtained at the Manna 
Kea and Cerro Paranal observatories. 

Finally, to obtain a good time coverage and increase the number of ob¬ 
servations we asked other observers to provide us with their images so that 
we could remeasure the astrometric positions prior to the encounter. We 
obtained and remeasured more than 200 images obtained at the following 
observatories: 

• Swift space telescope (this was hrst time Swift has been used to obtain 
small body astrometry); 

• RAS Observatory, Moorook; 

• Sutherland-LCOGT A and B; 

• iTelescope Observatory, Siding Spring; 

• Siding Spring-LCOGT A; 

• Blue Mountains Observatory, Leura; 
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• CAO, San Pedro de Atacama; 


Murrumbateman. 


For these stations, we performed the astrometry with track and stack pro¬ 
cedures, and corrected the astrometric positions by as much as 2". For each 
image, we estimated the astrometric error from the image quality, following 
a scale that assigned 1 pixel astrometric error to all the images where an 
obvious star-like central condensation was visible, and 2 pixels to more prob¬ 
lematic cases where the central condensation was either diffuse or blurred 
(due to poor seeing conditions and/or tracking). We evaluated a few in¬ 
termediate cases at 1.5 pixels. We set the data weights by conservatively 
applying a scale factor 2 to the estimated uncertainties. 

Figure shows the right ascension (RA) and declination (DEC)|^ residu¬ 
als against solution 102, which is our current best estimate for the comet’s 


trajectory and is described in Sec. |3.4[ It is clear how important the remea¬ 
surement and observation selection process was: the selected and remeasured 
astrometric observations have smaller errors than the excluded data. In par¬ 
ticular, some of the observations we hltered out during our selection process 
have residuals as large as 5", while the astrometric positions replaced by our 
remeasurements have residuals as large as 2". 


^In this paper RA and DEC are always referred to the equatorial J2000 reference frame. 
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2.2. Dynamical models 


To model nongravitational perturbations we first considered the classical 


formulation by Marsden et al. (1973): 


^NG = 9{.r){Air + Aat + Agh) 





( 1 ) 


where r is the heliocentric distance, m = 2.15, n = 5.093, k = 4.6142, tq = 
2.808 an, and a is such that g{l an) = 1 . Ai, ^ 2 , and are free parameters 
that give the nongravitational acceleration at 1 an in the radial-transverse- 
normal reference frame dehned byf = r/r, t = hxf, fi = rxv/|rxv|, 
where r and v are the heliocentric position and velocity of the comet. It is 
common practice to ignore the out-of-plane component, i.e., A^ = 0 . 

As of 2014 October 3, it was clear that nongravitational perturbations 
were needed to fit the observed data. We computed two different orbital 
solutions: 


• Solution 95, where Ai, A 2 , and A 3 are all determined as part of the 
least squares orbital fit; 

• Solution 97, where Ai and A 2 are determined from the orbital £t and 
A 3 is set to zero. 

Table [ 1 ] shows the values of the nongravitational parameters for solutions 95 
and 97, and Fig. shows the corresponding 6 -plane predictions. Solution 
95 corresponds to the more conservative approach where no assumptions are 
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Solution 


95 

97 

Al 

0 

1 

CD 

au/d^] 

18.8 ±4.8 

15.8 ±3.5 

A2 

0 

1 

CD 

au/d^] 

-11.2 ± 13.0 

4.4 ±3.5 

A3 

0 

1 

CD 

au/d^] 

-8.4 ±6.5 

0±0 


Table 1: Estimate of nongravitational parameters (along with Icr formal uncertainties) as 
of 2014 Oct 3. Solution 95 has all three nongravitational parameters Ai, A 2 , and as 
free parameters in the fit, while solution 97 is constrained to A 3 = 0. 


made on the nongravitational parameters. On the other hand, solution 97 is 
more aggressive, as is set to zero. The nominal value of for solution 


95 is larger than expected, e.g., the reference scenario in [Farnocchia et^ 
(2014) corresponded to a ±2 x 10“® au/d^ 3cr interval for A 3 . We thought 
that such a large value of A 3 was unlikely and could be due to still unresolved 
issues in the astrometry. Moreover, the A 3 detection was only marginal (l.Scr 
signihcance). Therefore, we decided to set A 3 to zero and choose solution 97 
as baseline ephemeris for the 2014 October 7 HiRISE observations described 
in Sec. 12.31 


2.3. HiRISE observations on 2014 October 7 

The uncertainties in nongravitational perturbations were too large to 
make sure the comet would be inside the HiRISE held of view for the close 
approach observations. For instance, if the pointing was centered on solution 
97 then solution 95 would be outside the held of view at closest approach. 

To tackle this problem, HiRISE observations from the Mars Reconnais¬ 
sance Orbiter took place on 2014 October 7. The much smaller range distance 
between C/2013 Al and MRO as well as the diherent geometry provided 
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much more leverage than Earth-based observations. The date for these ob¬ 
servations was chosen so that 1) the plane-of-sky uncertainty would be small 
enough to ensure that C/2013 Al would be in the held of view; 2) there 
would be enough time to transmit the images to the ground, perform the 
astrometric reduction, update the ephemeris and the pointing for the close 
approach observations; 3) the observations would provide the needed leverage 
to constrain the trajectory of C/2013 Al. 

The HiRISE observations utilized the three central broadband RED CCDs 


in the instrument focal plane (McEwen et ah, 2007), yielding images of the 
star held surrounding C/2013 Al that were ~20' on a side as shown in 
Fig. The scanning rate was selected to obtain the maximum useful inte¬ 
gration time of the HiRISE instrument, 2.56 seconds. The individual CCD 
images were corrected for dark current, instrument noise and for low fre¬ 
quency spacecraft jitter. The CCD images were joined together utilizing 
reconstructed SPICE kernels provided by NASA’s Navigation and Ancillary 


Information Facility (Acton, 1996) and SPICE related utilities in the Inte¬ 


grated Software for Imagers and Spectrometers software package (Torson and 


Becker, 1997), resulting in a single mosaic of the three CCDs to which an 


astrometric solution could be applied. To determine the background star 
positions on the CCD mosaic, their centroids were measured and a plate so¬ 
lution computed from a least squares ht of these star positions. The position 
of the comet was determined by measuring the centroid of the comet image 
and converted into RA and DEC via the plate solution. Since HiRISE is 
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a scanning instrument, the time of the observation was computed from the 
image line determined from the centroid measurement. Figure shows the 
orientation and scan direction of the three HiRISE images acquired on 2014 
October 7. 

The measured astrometric positions are reported in the hrst block of 
Table To account for the uncertainties resulting from centroiding, MRO 
ephemeris errors, and the small number of reference stars in the held of 
view, we conservatively weighted these data at 1". We chose these weights 
according to an assessment of the internal consistency of the data and the 
estimated systematic errors. 

Figure shows the HiRISE astrometric positions with respect to solu¬ 
tions 95 and 97. These positions are clearly inconsistent with solution 97, 
which we had deemed as more likely. On the other hand, solution 95 pro¬ 
vide an excellent prediction for the HiRISE observations, thus conhrming the 
(unexpectedly) large out-of-plane component of the nongravitational pertur¬ 
bations. 

The 2014 October 7 HiRISE observations removed the prejudice on 
by showing that the large out-of-plane perturbation seen in solution 95 was 
real and not an artifact of unresolved issues in the astrometry. Since solution 
95 and 97 would have led to different held of views for the close approach ob¬ 
serving sequence, these HiRISE observations were decisive to obtain accurate 
ephemerides for the Mars encounter. 
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Figure 3: The orientation and scan directions of the three 2014 October 7 observations 
from HiRISE. The green, cyan, and red colored boxes are the footprints of the three 
observations plotted on top of a digital sky survey field. The colored arrows indicate 
the direction of the scan of the corresponding footprint. The approximate position of 
C/2013 Al is indicated in white. 
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Date 


RA 

DEC 

ARA 

ADEC 

UTC 

[hh mm ss] 

[dd 

mm ss] 

n 

n 

2014-10-07.6273801505 

02 44 18.0206 

-15 25 39.461 

+0.212 

+0.100 

2014-10-07.6447532755 

02 44 18.1546 

-15 25 27.142 

+0.006 

-0.028 

2014-10-07.6621002778 

02 44 17.2358 

-15 25 16.345 

-0.630 

-0.099 

2014-10-17.2864254282 

02 

44 

16.40 

-14 

54 

09.3 

-0.038 

-0.049 

2014-10-17.3072505093 

02 

44 

15.54 

-14 

52 

34.9 

+0.084 

-0.046 

2014-10-17.4421608565 

02 

44 

16.16 

-14 

51 

30.4 

+0.093 

+0.670 

2014-10-17.4629921412 

02 

44 

15.23 

-14 

49 

48.3 

+0.443 

+0.422 

2014-10-17.5979258218 

02 

44 

15.84 

-14 

48 

29.8 

-0.372 

+0.523 

2014-10-17.6187525000 

02 

44 

14.83 

-14 

46 

38.4 

-0.265 

+0.614 

2014-10-17.7537702431 

02 

44 

15.54 

-14 

45 

01.1 

+0.024 

+0.284 

2014-10-17.7746015856 

02 

44 

14.40 

-14 

42 

58.6 

-0.113 

+0.298 

2014-10-17.9094922338 

02 

44 

15.16 

-14 

40 

57.8 

+0.128 

-0.049 

2014-10-17.9303112616 

02 

44 

13.88 

-14 

38 

41.7 

-0.444 

+0.368 

2014-10-18.0652698032 

02 

44 

14.69 

-14 

36 

09.1 

-0.241 

+0.252 

2014-10-18.0860832407 

02 

44 

13.32 

-14 

33 

37.2 

-0.047 

+0.394 

2014-10-18.2210950694 

02 

44 

14.15 

-14 

30 

22.5 

-0.072 

+0.451 

2014-10-18.2418975810 

02 

44 

12.60 

-14 

27 

30.8 

-0.246 

+0.126 

2014-10-18.3768168866 

02 

44 

13.49 

-14 

23 

19.4 

+0.336 

+0.075 

2014-10-18.3976231829 

02 

44 

11.73 

-14 

20 

01.1 

-0.211 

+0.084 

2014-10-18.5325476042 

02 

44 

12.65 

-14 

14 

28.8 

+0.190 

+0.340 

2014-10-18.5533507060 

02 

44 

10.71 

-14 

10 

35.9 

+0.721 

+0.435 

2014-10-18.6944945486 

02 

44 

12.99 

-14 

01 

59.8 

+0.472 

+0.138 

2014-10-18.7152678472 

02 

44 

05.55 

-13 

57 

05.5 

-0.456 

+0.032 

2014-10-18.8507619560 

02 

44 

11.78 

-13 

46 

21.1 

+0.109 

+0.122 

2014-10-18.8715381481 

02 

44 

02.65 

-13 

40 

15.2 

-0.236 

-0.495 

2014-10-19.0070050579 

02 

44 

10.13 

-13 

24 

17.5 

+0.262 

+0.002 

2014-10-19.0277652778 

02 

43 

58.53 

-13 

16 

19.1 

-0.016 

+0.029 

2014-10-19.1598250347 

02 

44 

07.04 

-12 

52 

32.7 

+0.085 

+0.517 

2014-10-19.1805583565 

02 

43 

56.30 

-12 

41 

15.2 

+0.059 

+0.343 

2014-10-19.3161075926 

02 

44 

02.73 

-11 

57 

00.0 

-0.134 

-0.253 

2014-10-19.3368152315 

02 

43 

47.09 

-11 

38 

56.9 

-0.076 

-0.267 

2014-10-19.4724035764 

02 

43 

54.07 

-10 

03 

09.3 

+0.114 

-0.041 

2014-10-20.2152129167 

14 

45 

14.70 

+19 

18 

44.9 

+0.324 

-0.639 

2014-10-20.3715103356 

14 

45 

00.17 

+18 

20 

00.8 

+0.255 

-0.683 


Table 2: C/2013 A1 astrometry obtained by the HiRISE camera of the Mars Reconnais¬ 
sance Orbiter. The last two columns contain the observed — computed residuals with 
respect to solution 102, described in Sec. |3.4^ 3 The RA residuals account for the spherical 
metric factor cos(DEC). 







Figure 4: Plane of sky on 2014 October 7 as seen from the Mars Reconnaissance Orbiter. 
The origin is set on the prediction corresponding to solution 97. Ellipses correspond to 
the 3cr uncertainty level. Dashed curves represent the plane-of-sky motion: C/2013 A1 is 
moving North-East. 
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Epoch TDB 

2014 Feb 3.0 

Eccentricity 

1.0007424 ± 0.0000029 

Perihelion distance 

1.3986359 ± 0.0000024 an 

Time of perihelion TDB 

2014 Oct 25.32823 ± 0.00019 d 

Longitude of node 

300.976279° ± 0.000048° 

Argument of perihelion 

2.44097° ± 0.00015° 

Inclination 

129.027341 ± 0.000020° 

Al 

(20.8 ± 2.6) X 10-y au/d^ 

A 2 

(-13.2 ±2.2) X 10-9 au/d^ 

A 3 

(-7.4 ± 1.3) X 10-9 au/d^ 

Mars encounter distance 

140751 ± 175 km 

Epoch of closest approach to Mars TDB 

2014 Oct 19 18:28:38.5 ± 1.2 s 


Table 3: JPL solution 101. The orbital elements are ecliptic heliocentric and error bars 
correspond to la formal uncertainties. 


3-4 ■ Ephemerides 

By including the 2014 October 7 HiRISE astrometry as well as some ad¬ 
ditional ground-based astrometry through 2014 October 13, which was the 
deadline for the hnal ephemeris delivery, we computed JPL solution 101. Ta¬ 
ble 1 ^ shows the corresponding osculating orbital elements, nongravitational 
parameters, and close approach distance and time. Figure shows the cor¬ 
responding 6 -plane prediction. 

Solution 101 was used to successfully observe C/2013 Al during the close 
approach to Mars. The data analysis of the acquired image^is still ongoing, 
though some preliminary results were reported by|Delamere et ah (2014). 


^Some of the images are available at http://mars.nasa.gov/comets/sidingspring/images 
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3. Post-encounter reconstruction 


3.1. Astrometry during and after the encounter 

Besides valuable information on the physical properties of C/2013 Al, the 
HiRISE observations obtained during the close approach to Mars provide fur¬ 
ther constraints to the comet’s trajectory. The series of observations taken 
at closest approach were assembled in a similar fashion to those of the 2014 
October 7 HiRISE observations. Scanning rates and integration times were 
varied to account for the relative motion of the comet and the spacecraft as 
well as the increasing brightness of the comet as the distance closed between 
the spacecraft and the comet. We could not obtain reliable astrometric mea¬ 
surements for a 20 hour period during closest approach as the scanning rates 
were too high to allow for accurate measurement of background stars. The 
second block of Table shows the measured astrometric positions. Similarly 
to the 2014 October 7 observations, these data were weighted at 1", with 
the exception of the last October 18 observation, which was weighted at 2" 
because of larger RMS errors of the plate solution. 

On 2014 October 19 and 20, C/2013 Al was observed by Hubble Space 


Telescope (HST, Li et ah, 2014b). We extracted accurate astrometry from 


drizzle-corrected (Hack et al., 2013) HST frames tracked on the comet. Since 
the reference stars in the held were trailed, we used trail-htting techniques to 
obtain proper astrometric solutions, referred to the PPMXL catalog ([Roeser 


et ah, 2010). In each image, the comet was clearly detected, and the central 


condensation was sharp and easy to identify, allowing high-precision astro- 
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metric measurements. To minimize the error from coma asymmetries, we 
used a small 2-pixel astrometric aperture in all cases. 


Because of the high quality of the HST images, the uncertainties of the 
astrometric positions were within 0.05". However, when htting these data 
we found residuals as large as 0.25". The cause of these larger residuals is 
the limited accuracy of the HST ephemeris. The geocentric coordinates of 
HST that we used are obtained from JPL’s Horizons system ( [Giorgini et al.| 
1996) and are derived from the Two-Line Element (TLE)|^ orbit solutions 


provided by the Joint Space Operation Centeij^ TLEs are not Keplerian 


elements and are hts to a particular dynamical model called SGP4 (Hoots 


and Roehrich| 1980). We analyzed the day-to-day variations by comparing 


the ephemerides of the HST TLE releases between 2014 October 2 and 17. 
The difference between the analyzed solutions suggest that the accuracy of 
the HST ephemeris is no better than tens of km. Therefore, we weighted the 
HST observations at 0.3". 

To extend the observed arc and include astrometry after the Mars flyby, 
we observed G/2013 Al from the Manna Kea observatory from February to 
April 2015. While the astrometric positions we reported to the Minor Planet 
Genter are from a 2.0 pixel (0.88") radius aperture, here we adopted the 
same linear extrapolation of the trend to a theoretical zero-sized aperture 


to compute the position of the cometary nucleus as desciribed in Sec. 2.1 


http://spaceflight.nasa.gOv/realdata/sightings/SSapplications/Post/JavaSSOP/SSOP_Help/tle_def.html 
®https://www.space-track.org/auth/login 
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Photocenters were computed using synthetic aperture radii ranging from 2.0 
to 6.0 pixels in steps of 0.1 pixels, so the extrapolation from 2.0 to 0.0 pixels 
represents 50 percent of the range actually £t. An average value for the 
tailward bias was computed from the trend seen in multiple exposures on each 
night (see Table |^. By default, the symmetric htting function computes the 
background level from an annulus surrounding the synthetic aperture with 
inner and outer radii of 2.5 and 5.0 times the radius of the synthetic aperture, 
respectively. For a pixel size of 0.439 arcsec, the background annulus was 
superposed on the coma of the comet, so some of the light from the comet 
was being subtracted from the area inside of the synthetic aperture. As such 
the reported photometry underestimates the true brightness of the comet 
for that aperture size. However, the goal here was accurate astrometry, not 
accurate photometry. We conservatively weighted these observations at 0.3" 
to 0.5", depending on the size of the bias corrections and the number of 


observations in the same observing night (Farnocchia et ah, 2015). 


3.2. Discrepancy between pre-encounter orbit and post-encounter data 

Though solution 101 provided successful ephemeris predictions for the 
HiRISE close approach observations of C/2013 Al, our attempt to extend the 
data arc revealed a clear discrepancy between the post-encounter data and 
the orbital predictions. For instance, when htting Ai, A 2 , and A 3 to the full 
data arc, we hnd A 3 = (—0.77±0.15) x 10“® au/d^, which is not statistically 
compatible with the much stronger negative acceleration measured before 
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Date 

East bias 

North bias 

2015-02-17 

0.175" 

0.025" 

2015-02-26 

0.175" 

0 .00" 

2015-02-27 

0.16" 

0 .00" 

2014-03-19 

0.145" 

0.04" 

2015-03-20 

0.155" 

0.035" 

2015-03-24 

0.13" 

0.045" 

2015-03-26 

0.145" 

0.055" 

2015-04-22 

0 .10" 

0 .12" 

2015-04-28 

0.06" 

0.125" 


Table 4: Astrometric corrections for the Mauna Kea post-encounter astrometry. These 
corrections are obtained by extrapolating to the limit of zero-sized aperture and were 
added to the positions reported to the MFC. 

the encounter (see Table . 

To quantify the discrepancy between the pre-encounter solution and the 
post-encounter data, we measured how much extending the data arc after 
the encounter changes the orbital solution. As shown by the fourth row 
in Table there is a statistically unacceptable >10cr correction from the 
pre-encounter best-£t solution to the full arc best-£t solution. 


3.3. Dynamical models 

The observed discrepancy points to the need of rehning the nongravita- 


tional perturbation model. The Marsden et ah (1973) model does a good job 
at htting the pre-encounter observations of C/2013 Al, actually none of the 
other models discussed below improves the £t to the pre-encounter arc. How¬ 
ever, this model proves inadequate as the trajectory is further constrained 
by post-encounter data. 
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We now analyze different nongravitational perturbation models, whose 
results are summarized in Table Besides the orbital £t statistics, to quan¬ 
tify the predictive power of the models Table shows the a shift applied to 
the pre-encounter solution to £t the post-encounter data. 


As described by Yeomans and Chodas (1989), the hrst rehnement is to 
introduce an offset AT of the activity relative to perihelion, which is where 


the Marsden et ah (1973) nongravitational perturbations peak. This asym¬ 
metry can be obtained by replacing g{r{t)) with g{r{t + AT)) in Eq. ([^. For 
C/2013 A1 the inclusion of AT does not result in any signihcant improve¬ 
ment: is essentially unchanged, the AT estimate from the orbital £t is 

1.6 ± 2.8 d, which is consistent with AT = 0, and we still have a lOa shift 
from pre-encounter best-£t solution to post-encounter best-£t solution. 

A higher fidelity model is the Rotating Jet Model (RJM), which is fully 


described in Chesley and Yeomans (2005) and we briefly recall here. The 
RJM considers the action of discrete jets on a rotating nucleus. Given a 
jet, we compute the mean acceleration averaged over a comet rotation. The 
acceleration depends on the orientation of the spin pole, e.g., dehned by RA 
and DEC, the jet’s thrust angle g, i.e., the colatitude of the location of the 
jet, and the strength Aj of the jet. The model can be extended to include a 
diurnal lag angle 6 and thus account for diurnal lag in the jet’s activity. All 
these parameters can be estimated as part of the orbital £t to the available 
astrometry. Note that the RJM does not depend on the rotation period and, 
unless 6 ^ 0°, is insensitive to the sense of rotation. 


24 











It is worth pointing out that Li et ah (2014a) detected two jet-like features 
in HST images obtained on 2013 October 29, 2014 January 21, and 2014 
March 11, which therefore support the application of the RJM. 

The RJM with one jet does not work well, e.g., is worse than that 
obtained with the Marsden et ah ( |1973 ) model. This behavior is somehow 
expected: as the comet moves around the Sun the latitude of the sub-Solar 
point changes and at some point the jet gets conhned in the comet’s shadow 
for the full rotation period, thus turning off the corresponding nongravita- 
tional accelerations. 

A possible rehnement to the single jet model is obtained by adding Ai, so 
that Ai gives the main nongravitational perturbation in the radial direction 
while the jet mostly contributes to the other components. This conhguration 
results in a signihcantly lower value and the a shift is 3.2, which is perfectly 
reasonable in an 11-dimensional parameter space (there is a 51% probability 
that Ax^ > 3.2^ with eleven degrees of freedom). 

The improvement obtained by adding Ai warrants the addition of a sec¬ 
ond jet to the model. After all, using two jets is a sort of generalization of 
the one jet -|- Ai solution. The orbital £t statistics and the a shifts of these 
two models are essentially equivalent. 

Further adding Ai to the two jets solution does not provide any improve¬ 
ment, thus suggesting that a third jet is not needed. On the other hand, the 
inclusion of the diurnal lag angle 6 lowers by 13 with respect to the best-£t 
two jets solution, which is a significant change for the addition of a single 
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parameter. Moreover, the two jets + 6 solution has the lowest a shift when 
the arc is extended after the encounter, which means that pre-encounter 
predictions are consistent with post-encounter data. This consistency can 
also be seen in the out-of-plane nongravitational perturbations, as shown 
in Fig. The large difference between the pre-encounter and the full arc 


best-£t solutions obtained with the Marsden et ah (1973) model is evident: 
the pre-encounter solution has much stronger out-of-plane perturbations, as 


already discussed in Sec. |3.2[ The RJM with two jets -|- 9 allows short-term 
variations due to the rapidly changing orientation of C/2013 Al with respect 
to the Sun. In particular, around the encounter one of the jets ends up in the 
comet’s shadow and becomes inactive. Thanks to the change of sign of the 
RJM out-of-plane perturbation, we can match both the large out-of-plane 
perturbations prior to the encounter as well as the much lower (in absolute 
value) average value over the full observed arc. 


In conclusion, the £t to the observational data and the HST images in Li 


et al. (2014a) warrant the usage of the RJM with two jets. Moreover, the fit 


also points to a diurnal lag in the jets’ activity. Interestingly, there is obser¬ 
vational evidence supporting the presence of such a lag for comets: measure¬ 
ments by the MIRO instrument on Rosetta show that for 67P/Churyumov- 
Gerasimenko there is a displacement between the sub-surface temperature 


peak and the sub-solar point (Gulkis et al., 2015). 
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Figure 5: Out-of-plane nongravitational perturbation for three orbital solutions: 1) Mars- 


den et al. (1973) model, best fit to the full arc, 2) Marsden et al. (19731 model, best fit to 


the pre-encounter arc, and 3) RJM model with two jets and diurnal lag angle, best ht to 
the full arc. 


3.4- Ephemerides 

As our best estimate of C/2013 Al’s trajectory we computed JPL solution 
102 by using the two jets -|- 6 model for nongravitational perturbations. 
Table shows the corresponding orbital elements, RJM parameters, and 
close approach data along with their formal uncertainties. 

Solution 102 provides a good ht to the astrometric data, e.g., the square 
root of the reduced is 0.29. In particular. Table contains the HiRISE 
astrometry residuals with respect to solution 102, which are well within the 
uncertainties used to weight the data, 2" for last October 18 observation, 1" 
otherwise. The left column of Fig. shows the RA and DEC residuals again 
solution 102 as a function of the observation epoch. 

Table [^reports the 6-plane coordinates ^ and ( for solution 102, which are 
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Epoch TDB 

Eccentricity 

Perihelion distance 

Time of perihelion TDB 

Longitude of node 

Argument of perihelion 

Inclination 

2014 Jul 18.0 

1.0008169 ±0.0000016 

1.39866662 ± 0.00000071 an 
2014 Oct 25.318385 ± 0.000096 d 
300.977577° ± 0.000026° 
2.438202° ± 0.000060° 
129.027337 ±0.000031° 

RA of spin pole 

62.8° ±7.0° 

DEC of spin pole 

14° ± 14° 

Strength of jet 1 

(86.5 ± 5.2) X 10-9 au/d^ 

Strength of jet 2 

(153 ± 67) X 10-9 au/d^ 

Thrust angle of jet 1 

60.2° ±3.1° 

Thrust angle of jet 2 

152.4° ±4.8° 

Diurnal lag angle 

9.4° ± 1.6° 

Mars encounter distance 

140496.6 ± 4.0 km 

Epoch of closest approach to Mars TDB 

2014 Oct 19 18:28:33.99 ±0.85 s 


Table 6: JPL solution 102. The orbital elements are ecliptic heliocentric and error bars 
correspond to la formal uncertainties. 



also shown in Fig. as well as those obtained with other dynamical models. 


The formal nncertainties in ^ and ( are ~4 km, while using the Marsden et 


al. (1973) nongravitational perturbation model can cause differences in as 


large as ~10 km. Despite the small encounter distance, 140500 km, the Mars 
flyby caused modest changes to the orbit of C/2013 Al because of the high 
encounter velocity of 56 km/s and the relatively small mass of the planet. 
For instance, the angle between the incoming and outgoing relative velocities 
of the comet was only 0 . 01 °. 

Interestingly, solution 102 comes with an estimate of C/2013 Al’s spin 
pole. While this pole allows us to £t the data and accurately describe the 
comet’s trajectory, it is important to point out that this estimate is model 
dependent and strongly relies on the accuracy of the RJM and its assump¬ 
tions, e.g., number of jets, simple rotation, etc. Table shows the different 
estimates of the spin pole with the Marsden et al. ( |1973 ) model (spin pole 
approximated by the direction of the (Ai, A 2 , A 3 ) vector at time AT, as done 


m 


Chesley and Yeomans, 2005) and the different conhgurations of the RJM. 


In particular, setting the diurnal lag angle to 0° moves the pole by 36°. 

Comet C/2013 Al is still observable in the northern sky. Figure |^shows 
the predicted total magnitude and declination until 2020, when C/2013 Al 
will be at more than 15 an from the Sun. The next few years are the 
last chance of observing the comet for quite some time. In fact, the post- 
perihelion orbit of C/2013 Al with respect to the Solar System barycenter 
is elliptic and has an orbital period of more than 500000 years. 
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Figure 6: Predicted magnitude (solid line) and declination (dashed line) as observed from 
the geocenter. Grayed parts correspond to solar elongation < 60°. 

4. Conclusions 

The campaign for observing C/2013 Al from the Mars orbiting space¬ 
crafts during the 2014 October 19 encounter was challenging in terms of 
computing reliable comet trajectories. The small held of view of Mars Re¬ 
connaissance Orbiter’s HiRISE camera and the close encounter distance set 
tight requirements on the ephemeris accuracy. 

One of the hrst problems we had to face was the generally poor quality of 
C/2013 Al’s astrometric observations. In fact, comet astrometry is often af¬ 
fected by tailward biases as the center of light does not generally correspond 
to the brightness peak. To mitigate this issue, we hrst applied a strict hlter 
to the astrometric data by only selecting observations coming from profes- 
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sional surveys. Then, we added observations for which we directly performed 
the astrometric reduction either because we directly performed the observa¬ 
tions or the original observers provided us with their astrometric images. 
This process allowed a more stable orbit update sequence as new data were 
available. 

Nongravitational accelerations often are the limiting factor in terms of 
comet ephemeris accuracy. To constrain the nongravitational perturbations 
acting on C/2013 Al we planned observations from the Mars Reconnaissance 
Orbiter HiRISE camera twelve days before the Mars encounter. Thanks to 
the favorable geometry, these observations were decisive to secure the trajec¬ 
tory of C/2013 Al and to provide accurate ephemerides for the encounter. 
In particular, the HiRISE observations clearly revealed a signihcant out-of¬ 
plane nongravitational acceleration thus showing how the common practice 
of restricting to in-plane perturbations can lead to inaccurate predictions. 

Though we only have astrometric data for a single revolution, the good 
observational coverage and the strong leverage provided by the HiRISE obser¬ 
vations acquired as the comet flew by Mars showed that the classical [Marsdeii 


et al. (1973) model for nongravitational perturbations is not accurate enough 


to model the trajectory of C/2013 Al over the full observed arc. We there¬ 
fore resorted to the Rotating Jet Model (RJM), which is a higher hdelity 
model that allows us to nicely £t the full observed arc thanks to the action 
of two jets. Moreover, the RJM provides an estimate of the spin pole of 
C/2013 Al. Though the reliability of this estimate strongly depends on the 
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accuracy of the model, Chesley and Yeomans (2005) already used the RJM 


to quite successfully estimate the spin pole of 67P/Churyumov-Gerasimenko: 


the difference between prediction and the actual spin pole orientation (Sierks 


et ah, 2015) is only 13°. The analysis of the images acquired during the close 


encounter will possibly provide an independent estimate of the spin pole to 
be compared to that obtained from the orbital dynamics. 

Future work will be related to further rehnements of the nongravitational 


perturbation. For instance, the standard g{r) function used by the Marsden 


et ah (1973) model and the RJM is driven by the sublimation rate of water ice 


as a function of heliocentric distance. Directly estimating the H2O and CO2 
production rates of C/2013 Al and future astrometric observations could 
allow one to derive a more C/2013 Al specihc g{r). 
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